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Abstract 

We consider the one-dimensional spinless Falicov-Kimball model of itinerant fermionic 
particles ("spinless electrons"), which can hop between nearest-neighbour sites only, and 
of immobile particles ("classical ions"), with an on-site attraction. Extensive studies of 
the ground-state phase diagram of this system and its higher dimensional counterparts, 
carried out up to now, concentrated on determining ground-state arrangements of ions 
on the underlying lattice, while the properties of electrons were typically ignored. We 
report studies of short- and long-range correlations between electrons, and between ions 
and electrons, and of the spatial decay of electron correlations (decay of single-particle 
density matrix), in the ground state. The studies have been carried out analytically 
and by means of well-controlled numerical procedures. In the case of period 2 ground 
state, the single-particle density matrix has been expressed in terms of a hyper geometric 
function, and its spatial decay has been extracted. Numerical calculations have been 
done for open chains of various lengths (up to a few thousand sites), in order to control 
the chain-size dependence of correlations and to extrapolate the results to the limit of 
infinite chain. A part of the obtained results refers to tight-binding electrons subjected 
to a periodic external potential due to the ions, which constitute simple models of metals 
and insulators. 
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1 Introduction 



The Falicov-Kimball model emerged from solid-state theory in 1969, as a simple model for 
semiconductor- metal transitions [jl|. Since that time it has become an important standard 
tool for studying various properties and phenomena that occur in strongly correlated fermionic 
systems, ranging from metal-insulator and mixed- valence phenomena to the Peierls instabilities 
and crystallization. For an overview and an extensive list of relevant references the reader is 
advised to consult 

In the present paper we are studying the simplest, skeleton version of the model proposed 
by Falicov and Kimball, which is usually referred to as the spinless Falicov-Kimball model. This 
is a lattice system that consists of two sorts of spinless fermions, called electrons and ions. The 
electrons can hop between nearest-neighbour sites only, while the ions are immobile. The only 
interaction occurs between the electrons and ions and it is an on-site attraction of strength U. 
Throughout this paper, the underlying lattice is one- dimensional (later on called the chain), 
so it can be identified with integers Z. When the chain is finite and consists of L sites labeled 
X = 0, . . . , L — 1, we impose boundary conditions b, which in this paper are either free {b = f) 
or periodic {b = p). Then, the whole two-component system is governed by the Hamiltonian: 

L-l 

Hb = - ^ (tb,xa^a:^+i + tb,xa^^iax + Uw^a^a^) , (1) 

x=0 

where the creation and annihilation operators of an electron at the site x of the 

chain, tf,^x is the hopping rate between the sites x and x + 1 (which will be set site-independent, 
except at the boundary of a chain), and Wx stands for the occupation-number operator of an 
ion at site x. Since the occupation-number operators commute with the Hamiltonian (|lD, 
we can identify them with classical variables taking values zero or one. Consequently the total 
number of ions, Ni = ^^-Wx? is conserved. Moreover, the total particle-number operator of 
electrons, A'^g = Ylx^t^x, commutes with Hb, so the total number of electrons, Ne, is also 
conserved. Therefore, we can study the equilibrium properties of our system in the canonical 
ensemble, with a fixed number of sites L, the electron number Nf, and the ion number iV,. The 
corresponding canonical partition function reads: 

Zb= Yl Tr^eexp(-/5i/,), (2) 

where w stands for the function that to a site x of the underlying chain assigns the value Wx = 
0, 1, which can be viewed as a sequence w = WqWi . . . (hereafter called the ion configuration) , 
the trace is over the Hilbert space of A^^e electrons, and (3 is the inverse temperature. 

Despite the apparent simplicity of the form of Hamiltonian (|l]) and its counterparts on 
higher- dimensional lattices Z"', d > 1, the question whether the system is capable of describing 
cooperative effects had remained open for almost two decades. It was affirmatively answered in 
1986 by demonstrating that the half-filled system exhibits a staggered type of long-range 

order in the ionic subsystem at zero temperature, and showing that this staggered order persists 
at sufficiently low temperatures in two- and higher- dimensional systems 0. 

Since then, one can observe an increased interest in the model and many properties of 
the Falicov-Kimball model have been established by means of rigorous methods as well as by 
means of various sorts of numerical procedures. Not surprisingly, a vast majority of obtained 
results refers to the ground state (see 0, |^, ^, and references quoted there). As there are 
good reasons to believe that studies of the ground state of the one-dimensional model may 
reveal interesting cooperative phenomena, which herald the appearance of similar cooperative 
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phenomena in higher- dimensional systems at zero and sufficiently low temperatures, a lot of 
efforts at determining ground-state properties have been made |Q]-|T^. 

As a matter of fact, all the attention, not only in the last quoted papers but also in many oth- 
ers whose references can be found in ^, ^ , has been concentrated on the ionic subsystem. 
Specifically, the most favourable energetically ionic configurations and ionic correlations, at 
different values of the electron-ion coupling parameter U, have been determined. The problem 
of correlations between the itinerant fermionic particles had not been usually addressed. To the 
best of our knowledge, there are just a few exceptions. As mentioned in [0] the results on decay 
of superconducting correlations, obtained in |]TU|, ^ for general itinerant-electron systems, can 
be applied to the ground state of the one-dimensional Falicov-Kimball model (p. In particular, 
it follows that the ground-state average of a^ay decays exponentially with the distance between 
the sites x and y, i.e. there is no so called off-diagonal long-range order. Similar results, but 
for multi-point products of fermionic creation and annihilation operators (that cannot be ex- 
pressed in terms of occupation- number operators), have been obtained in the strong coupling 
regime in In p2| a class of models related with the spinless Falicov-Kimball model, with 



an emphasis on the static Holstein model, has been analyzed rigorously. As a byproduct of 
considerations carried out for the static Holstein model, the authors have obtained a relation 
between grand-canonical averages of one-site ion- and one-site electron-occupation numbers at 
half-ffiling. This result seems to be the first one referring to correlations between the electrons 
that can be expressed as an average of electron occupation-number operators. The mentioned 
relation implies that, at half-filling, the staggered long-range order in the ion subsystem is 
accompanied by the same kind of long-range order in the electron subsystem. 

It should be emphasized here that we are concerned only with low-dimensional systems. 
There is an extensive literature concerning the limit of infinite dimensions, where calculations 
of some correlation functions, like linear-response functions, are quite common (see for instance 

il). 

The purpose of this paper is to examine the properties of the electron subsystem of the 
model (|l|), at zero temperature. It is achieved by means of a number of short- and long-range 
correlation functions of electrons, and electrons and ions, that can be represented as canonical 
averages of operators constructed out of occupation-number operators of electrons and ions. 

The paper is organized as follows. We start, in Section 2, with a brief discussion of canonical 
ground-state averages. In Section 3, we define and discuss some general properties of short- 
and long-range correlation functions, to be studied in the subsequent sections. Then, in Section 
4, we present analytical results for closed finite chains and for infinite chains. After that, in 
Section 5, we give numerical results, obtained by means of numerical exact diagonalization. We 
end up, in Section 6, with a discussion and a summary of the obtained results. 



2 The canonical ground state and ground-state averages 

There are two approaches that enable one to study ground-state properties of a system. The 
first one is a quantum-mechanical approach, usually adopted in the papers concerning ground- 
state of Falicov-Kimball models, and consists in considering a finite (and then an infinite) 
system exclusively at zero temperature. The second one is a thermodynamic approach, where 
the considerations start with nonzero temperatures and then the limit of zero temperature is 
taken. 

In the quantum-mechanical approach, one makes use of the fact, mentioned in Introduction, 
that a system of electrons and Ni ions with Hamiltonian (|l]) , confined to a chain (of L sites 
and with boundary conditions b), can be considered for a fixed ion configuration w. The 
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corresponding Hamiltonian, denoted Hf,{w), describes tight-binding electrons in the external 
potential —Uwx- In this case, the ground state of Ng electrons, \w)b, that is the eigenvector of 
Hbiw) to the lowest energy Eb{w), is the so called Fermi-sea state. Then, the ground state of 
a finite chain with A^e electrons and Ni ions is given by the set Gb of finite-chain ground-state 
ion configurations gb and the corresponding Fermi-sea states \gb)b of energy Eb{gb), such that 

Eb{gb) = min Ebiw). (3) 

To get thermodynamically relevant quantities, we first introduce the finite-chain ground-state 
energy densities: 

eb{w) = L~^Eb{w), Cbigb) = L~^Eb{gb) = min eb{w), (4) 

for a given ion configuration w, and for given particle numbers N^,, N^, respectively. Then, 
we go to the thermodynamic limit, denoted lim^^oo; i-e. we send L to infinity in such a way 
that the both ends of the chain (if any) become remote from any fixed site, and L~^Ne — > pe 
- the electron density and L"^Ni ^ pi - the ion density, for some Pe,Pi G [0,1]. In the 
thermodynamic limit and for a fixed ion (infinite-chain) configuration w, the (infinite-chain) 
ground-state energy density is 

e{w) = lim eb{w), (5) 

L—*oo 

and is independent of boundary conditions. In the r.h.s of (|]), the ion configuration w stands 
for the restriction of the infinite-chain configuration w to the finite chain. For given particle 
densities pe, Pi, the infinite-chain ground-state energy density, e, is 

e = mine(w) = e{g), g E G, (6) 

where the minimum is taken over infinite-chain configurations of ions and the set G, defined 
by (|^), stands for the set of infinite-chain ground-state ion configurations. Alternatively, the 
ground-state energy density can be obtained as 

e = lim Cbigb), gt e Gb. (7) 

L—fOo 

In the thermodynamic approach, we start with nonzero temperatures and define the finite- 
chain internal-energy density cp^b, 

ei3,b = L~^{Hb)i3,b, (8) 

where 

{Hb)f3,b = Z-' Yl T^^N^expi-pHb) (9) 

is the canonical average of Hb. Then, the (infinite-chain) internal-energy density, e/3, is given 
by = limj;^^oo ^i3,b, and is independent of boundary conditions. Thermodynamically, the 
ground-state energy density, e, equals to e = lim/3^00 C/j. Reversing the order of the limits 
L ^ 00, /9 — >^ 00, we recover the formula (|^) of the quantum- mechanical approach 

e = lim Cbigb), eb{gb) = lim C/s^b = L'^b{gb\Hb\gb)b, gb e Gb. (10) 

L— ►oo /3— >oo 

To define a ground-state average, we start with nonzero temperatures. Let A be an operator 
built out of a;^, and/or w^, depending on a finite number of sites only. By definition, the 
finite-chain canonical average of A, {A)p^b, is 

{A)p,b = Z,' Yl Tr^A^M-PHb). (11) 



4 



The infinite-chain hmit of amounts to the canonical average of A, {A)p, 



hm {A)p^h = {A}^, 



(12) 



which is independent of boundary conditions and translation invariant. Then, the infinite- chain 
ground-state average of A, {A) , is defined as 



lim {A)p = {A). 



(13) 



The ground-state average (A) shares the above mentioned properties of (A) p. 

By rigorous results of |jl^ and numerical results of [1^, |13|, [l^ , for pairs of rational densities 
(pe,Pi), satisfying a linear relation, and for suitable values of electron-ion attraction U, the set 
G of ground-state configurations of ions consists of a finite number of elements. This happens, 
for instance, if pe = p/q, where the natural numbers p and q are relatively prime, and pi = pe 
(the so called neutral case), for sufficiently large U |I2|, |13|. Then, G consists of q periodic 
configurations of period q (the so called most homogeneous configurations) that differ from 
each other by a translation. In such a case, the ground-state average of A equals to 



(14) 



where g in \g)\, stands for the restriction of a (7 G G to the finite chain. As a finite-chain 
approximation to (A), converging to {A) as L ^ 00, we choose {A)h given by 



i,A), = q-^Y.^{9\^\3^^- 



(15) 



In what follows, we refer to the ground-state average (^) as the symmetric average, while the 
quantum average i,{g\A\g)i, is called the broken-symmetry average. Since the quantum aver- 
ages b{g\A\g)b appear frequently throughout the paper, we introduce the abbreviated notation 
{A)h^g = h{g\A\g)b. In general, a finite set of Q ground-state configurations G may consist of 
configurations related by translations and refiections, with Q > q. 



3 Short- and long-range correlation functions: defini- 
tions, general properties, and the method of calcula- 
tion 

Throughout the sections that follow, the set G stands for the set of periodic ground-state 
configurations of the infinite chain, for given pe = p/q, with p and q being relatively prime 
natural numbers, pi = mpe, with m = 1,2,..., and for suitable values of U. Later on, the set of 
U for which G is the set of the ground-state configurations will be called the stability interval of 
G, and denoted AUg- The set G consists of q ion configurations, with period q, which are known 
as the atomic (m = 1) or m-molecular [m > 2) most homogeneous ion configurations 0, |12], |T3[| . 



Following |T^, the set of m-molecule most homogeneous configurations corresponding to the 
electron density pe = p/q is denoted [p/q]m- The elements of [p/q]m can conveniently be 
described by specifying their restriction to the unit cell, whose sites are labeled x = 0, 1, . . . , g — 
1, i.e. as a sequence g = g^gi . . . gq-i. The positions of ions in the unit cell of the representative 
configuration of [p/q\m can be obtained as the solutions kj to the equations: 

pkj = j mod q, j = 0, 1, . . . , mp — 1. (16) 
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The remaining q — 1 configurations of \p/q\m are obtained by translating the defined above 
representative configuration. 

Moreover, the number of sites in the chain, L, is set to be a multiple of q, so that L~^Ne = pe, 
and L~^Ni — pi. For the ion subsystem the correlation functions are defined in Appendix 
A. Below, we define analogous correlation functions for the electron subsystem. The hst of 
correlation functions, studied in the sections that follow, opens a one-point function, the local 
electron density in the ion configuration g: 

{n^)h,g, geG. (17) 

It is related to Pe by 

i^-' J](n.)6,g = Pe. (18) 

X 

For periodic boundary conditions {b — p), {nx)p,g has period q: {nx)p,g — {nx+q)p,g, and 

q-l 

Pe = {nx)p = q~^ ^{nx)p,g = J^^'^'^+y'^P'S- (1^) 

sec? y=Q 

Then, we consider a number of two-point correlation functions. The first one is the density- 
density correlation in the ion configuration g, {nxny)b,g, which (by the Wick's theorem) can be 
expressed as 

{nxny)b,g = {nx)b,g{ny)b,g - \ {atay)b,g\^. (20) 

For the periodic boundary conditions, {nxny)p,g has period q, {nxny)p is translation invariant, 
and the symmetric average {{rix — pe) {ny — pe))p can be written as follows 

{{rix - pe) {uy - pe))p = q~^ '^{nxny)p,g - pi 

geG 

seG 9 

The symmetric average of the above kind is used to define the electron- electron short-range 
correlation function, Cb,x{l), 

^h,x{l) = {{rix - Pe) {nx+i - pe))b- (22) 
In an analogous manner, we define the electron-ion short-range correlation function, Sb,x{l)'- 

\x{l) = q~^ "^{{gx - mUx) {gx+i - mnx+i))b,g- (23) 

In the case of the periodic boundary conditions, the correlation functions, £5 ,,.(/) and Sb^x{l), are 
translation invariant (i.e. depend only on /). To define the long-range correlations corresponding 
to the above short-range correlations, we start with the order-parameter operators Ok, Og, and 

Og,g'- 

6fe = L-^J]e^'=^(n,-pe), (24) 

x=0 
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L-1 



Og = L ^^{Qx- Pi) {rix - Pe) , (25) 



x=0 

and 

L-l 



Og,9' = L {g^ - Pi) {g'^ - mn^) . (26) 



x=0 



Taking broken-symmetry averages of the order-parameter operators, we obtain the correspond- 
ing order parameters: 



L-l 



x=0 



L-l 



{Og)b,g> = L ^"^igx- Pi) {rix - Pe)b,g'- (28) 



1=0 

and 

L-l 



{Og,g')b,g' = L ^^{gx- Pi) {g'x - mn^)b,g'- (29) 

x=Q 

The symmetric averages of the order-parameter operators vanish: 

= (Og), = q-' Yl (Og,o')b,9' = 0. (30) 
g'eG 

By means of the order-parameter operators we define the long-range correlation functions Vb{k), 

n{k) = (0+4)6 

L-l 

= L-2 J2 e^'^^-yHin, - Pe) (% - Pe))b, (31) 

x,y=0 



C,^q-'J2{0'g)b 
geG 

L-l 

= Yl ^^3x - Pi) {gy - Pi)] {{Ux - Pe) {Uy - Pe))b, (32) 

3:,y=0 geG 



and 



S, ^ Y (Olg\, 

g,9'eG 

L-l 

Y Yl ^^3x - Pi) i9y - Pi)] Y^^d'x - mn^) {Oy - mny))t„g'. (33) 

x,y=0 geG g'eG 



In the case of the periodic boundary conditions, the long-range correlation functions Vp{k), C 



and Sp assume the form 



L-1 



P,(A;) = L-i5^e-%(0, (34) 



1=0 



L-1 



£p = L-i5^Ep(/)£p(Z), (35) 



1=0 



L-1 



5p = L-i5^Ep(/)5p(/), (36) 



1=0 



where Cp{l) = Cp^x{l), <Sp{l) = Sp^xil), and Ep{l) is defined in Appendix A. The function 
k — i> Vh{k) is usually referred to as the static structure factor p5[. Since, the following relations 
hold, 



L-1 
X, y=0 

L-1 

e^'^^-y^\{a^a,),J', (37) 

X, y=0 

L-1 

{Ol)b,9' = L'"^ Y - pi) i9y - pi) {rix - Pe)b,g'{ny - pe)b,g' 

x,y=0 

L-1 

-L~^ Y - Pi) {gy - P^) (38) 

X, y=0 



and 



L-1 

{Olg')b,9' = Y ^3x - pi) {gy - Pi) {g'x - mnx)b,g'{g'y - mny)i,^g, 

X, y=0 

L-1 

~L'^ Y i9x- Pi)igy~ Pi)\{atay)b,g'\'^, (39) 



X, y=0 



and since, for periodic boundary conditions, the averages {Og)p, q ^J2g'£G('^'g,g')p'9'^ 
{0~lOk)p,g are translation invariant, hence independent oi g E G, 



L-1 



Vp{k) ^ {dtd,)p=\{d,)pj'-L-'J2^"''\(^tcix^i)p,y\', (40) 



=0 



^ 2 

,9' 



Cp^q-'J2(Ol)p = Q-'J2{(Og)p, 

g&G g'£G 
L-1 

-L-^ ^9x - Pz) {gx+i - p^) Y K^x«x+i)p,g'P, (41) 

(=0 g'eG 
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and 

g,9'GG 9'GG g'GG 

L-1 

-^"^ E '^^'^ ~ ^3x+i - Pi) q'^ E I ^^t(^x+i)p,9'?- (42) 

«=0 g'GG 

For the one- dimensional systems of the kind considered here, and in the hmit L ^ oo, the 
second term in (0), (|4l|), (|4^ ) vanishes (see [|19|, and the sections that follow). Thus, 



in the thermodynamic limit we find simple relations between long-range correlations and the 
corresponding order parameters 

V{k) = hm Vp{k) = hm {d+dk)p = I lim (6,%,,^, (43) 

L— >oo L^oo L— ►oo 

2 



hm £, = g-i V ( lim(6,)p,,,) , (44) 



g'€G 



and 



5 = hm S, = q-' J2 (li"^ {Og,g,),,S. (45) 



g'GG 

The relations (^), (0), ( ^5] ) hold for other boundary conditions as well. If the quantities 
V{k), C are strictly positive, one says that the considered system exhibits the long-range order. 

In order to calculate the correlation functions defined above, one can use the fact that for 
a fixed ion configuration w the Hamiltonian of the system, Hb{w), is a second quantized form 
of the one-electron Hamiltonian h^iw). If we denote by {\x)}x=q,...,l-i the orthogonal basis of 
one-electron states, such that creates an electron in the state \x), then the matrix elements 
of hb{w) in the basis {|x)}a,-=o,...,L-i are defined by 

L-1 

Hh{w) = ^{x\hb{w)\y)atay, (46) 

X, j/=0 

i.e. explicitly, all the non-vanishing matrix elements are given by 

{x\hb{w)\x) = ~Uwx, {x\hb{w)\y) = ~tb,x if y = x±l. (47) 

Let {|f be the orthonormal basis built out of the eigenstates of hb{w) to the eigen- 

values A^, such that < A„/ if f < v'. Then, the unitary matrix U, with the following matrix 
elements Uxv' 

Kxv = {x\v), (48) 

diagonalizes the matrix of hb{w), 

^U^y{y\h{w)\x)Ux^, =K5vv'- (49) 

Moreover, the set of operators {6^, 6^}, w = fo, . . . , wl-i, defined by 

L-1 

K = '^{v\x)ax, (50) 

x=0 
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satisfies the canonical anticommutation relations, and 



L-l 



X, y=0 V 



(51) 



Since, 



x\v)by, 



(52) 



we can express the basic electron-correlation functions, {rix), (a+Oj^);,^^, and {nxny)^^^, in terms 
of the site-components {x\v) of the eigenvectors \v): 

{n,.)b,w = b{w\a^ax\w)b = \{x\v)\'^, 



V<Vp 



{a^ay)b,yj = b{w\a^ay\w)b = ^ {v\x){y\v), 



V<Vp 



V<-Vp 



{nxny)b,w = b{w\a^axa:l^ay\w)b = ^ \{x\v)\^ ^ 

V<Vp 



V<.Vp 



(53) 



where vp stands for the label of eigenvectors \v), such that there is exactly A^^e eigenvectors with 
V < vp- Consequently, all the correlation functions defined in this section can be expressed in 
terms of the site-components {x\v) of the eigenvectors of hb{w). For some low-period ion config- 
urations the eigenproblem can be solved exactly, while for any other ground-state configurations 
of interest, and arbitrary boundary conditions, one can resort to numerical exact- diagonaliz at ion 
procedures. 



4 Correlation functions — exact results 

In this section, we calculate exactly ground-state correlation functions of a finite chain with 
the periodic boundary conditions only (so without any risk of confusion the subscript p by the 
averages can be dropped), and then the infinite-chain limits, in two cases. 

First, for completeness, when pe G [0, 1] and pi = 0. For such a particle densities, G consist 
of a single ion configuration, called the empty configuration, for any U. In this configuration, 
= for all X, so we denote it g = 0. 

Second, when p^ = 1/2 and p, = 1/2, and G = [l/2]i consists of the two checkerboard 
configurations: ch^ = 10 and ch^ = 01. This holds also for any nonzero U, i.e. At/[i/2]i ~ (0' ^) 
i- 

When gx = for all x, all the non-vanishing matrix elements of hp{0) are given by 

{x\hp{0)\x + 1) = -tp^x, tp^x = l, X = 0,1, . . . ,L - 1, (54) 

Let k G {27tI/L : / = 0, 1, . . . , L — 1}, so that for any L, < k < 27r. The eigenvectors of hp{0), 
denoted \k), are specified by their site components, {x\k), 

{x\k) = L-^/^e'^^, (55) 
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while the corresponding eigenvalues, e^, read 

Ek = —2 cos k. (56) 
Introducing kp = n^N^ — 1)/L, we obtain 

(«X+Oo = y E exp(^A;/) = i[l + 2 cos(fc/)] = l ^^^^f ^J^f . (57) 

L ^-^ L ^-^ L sm 7r//L 

k<kp, 0<k<kp ^ ' ' 

fc>27r — kp 

Thus, in the infinite-chain limit 



, / ^ V ,9 1 sin^(7rpe0 /^^N 
hm \{ata,+i)o\' = (58) 



lim (n^n^.+z)o = pI 



1 sin^(7rpe/) 



L— »oo 



/2 



(59) 



The analogous expressions for continuous counterparts of our system can be found in p5[] . 
The case of open boundary conditions has been considered in detail in p6|]. Let us mention 
that the expression for the density- density correlation function {n^n^^i), obtained in |2^, differs 
from our result (^), since the limit L ^ oo is constructed in a different way, namely as the 
limit of half-infinite chain (only one end becomes remote from a fixed site). 

More interesting is, of course, the case of the checkerboard ground-state configurations 
g = ch^, ch^, which are related by the primitive translation of the chain. In this case, all the 
non- vanishing matrix elements of /ip(ch^) are given by 



{x\hp(ch^)\x + 1) 



(60) 



In the basis {|fco); l^o + ^r), |A;i), |A;i + vr), . . .}, k = 2ttI/L, I = 0, 1, ... , L/2 — 1, obtained by 
reordering the basis the matrix of /ip(ch^) becomes block-diagonal. The diagonal blocks 

are 2 by 2 matrices [/ip(ch^)]^ that in the basis {\k), \k + vr)} assume the form 



— ^ 

U 



Thus, we obtain easily the eigenvalues, A±{k), 



U 



A^{k) = --±A, 



A 



_u 
-£k - 2 



si + a' 



U 
2 ' 



± 

chi 



and the corresponding eigenvectors, given by: 

[/i,(ch^)]Jfc)± =A±(fc)|fc) 
and 

\k)^K = ilkr' + Pk\k + ^)) , = ekTA, 



(61) 



(62) 



(63) 



where the normalizing factor 7^ is 



(7, 



±\2 



±\2 



(64) 
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Since in what follows, we shall be interested in pe < 1/2 only, we can restrict our considerations 
to the lowest L/2 eigenvalues A~(fc) and the corresponding eigenvectors \k)~ , and set 

\k)ch = \k)-f,, /3k= f3k = ek + lk = lk- (65) 
Then, the site-components of eigenvectors \k)ch read 

{x\k),h = 4^ - [ae^^^ + /?,e^('=+-)^) , (66) 
VL Ik 

and consequently, setting (ch^ja+a^ilch^) = {a^ax)i, the local density {nx)i reads 

K)i = (aX)i = 7 Yl \W)ch\'=Pe + {-iraTo{0), (67) 



where 



In the thermodynamic limit. 



1 f 

lim ro(/) = — dk 

L^oo 71 J 



cos(fc/) 



A 







■Kpe npe 

1 /" cos( kl) K /" cos( kl) 1 9 1 , 

' dk , ^ ^ = — / dfc , ^ ^ , Pe < -, = 777TTT < 1. 69 



Since, the function limL^ooTo(0 appears in many expressions for correlation functions, in the 
sequel, it is interesting to see its dependence on pe and U. This is shown in Fig. 0. The analytic 
expressions for the asymptotic behaviour of this function, for f/ — > and for U oo, are 
given in Appendix B. Using these asymptotic formulae we find the corresponding asymptotic 
expressions for the local density {nx)i: for U —>■ 

L^oo^ ' \ Pe + {-!) a{2n) ^ (ao + 21n2 - Ina) , Pe = ^, 
while for U ^ oo and any pe < 1/2, 

lim lim (n^)i = Pe + (-l)Ve, (71) 

U — >oo L— >oo 

which in the case of half-filling (i.e. for pe = 1/2) amounts to ch^. Moreover, in the case 
of half-filling, the local electron-density of an infinite chain can be expressed by the complete 
elliptic integral of the first kind K [^], at the point k^, 

tt/2 

, , 1 , ^^.Uk f dk 1 , ^^Uk , r.^ 1 , ^ 

- 2 ^ (-1)^ 4^ / TTTOT - 2 - (-'r^'<(^'^, P. ^ (72) 
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Now, consider the two-point correlation function {a'^ax+i)i = {ch.^\a'^ax+i\ch.^)ch- For < 
1/2, we find 

{a^a^+i)i= ^ i{k\x){x + l\k)i 

k>n—kp 

= z E Cw+{-iypi+{-irc.Pk[i+{-iy]}, (73) 

fc>7r— fc^ 

which for nonzero and even / assumes the form 

fc>7r — fcjjT 

= (4a.+^)o + i-lTa^ (71^ + 2 5^ ^ (a+a.+,)o + (-l)^«ro(2s), (74) 

while for odd the form: 

(«X+/)i = ;^ E =^1^^)' ^ = 25 + 1, 5 = 0,1,..., (75) 



where 

+ ^ E oos(«)(-|)}. (76) 

V -1 -r u; o<fe<feF J 

We note that the functions Tq and Ti satisfy the relation: 

ri(2s + l) = ro(2s + 2)+ro(2s). (77) 

Summarizing, for pe < 1/2, the two-point function {a^ax+i)\ reads 

. , _ / ri(0, / = 2s + l, s = 0,l,..., . . 

^ ^ I (aX+^)o + (-ir«ro(/), ; = 2s, s = l,2,.... ^'^^ 

As expected, since for f/ — > 0: Ti(/) ^ (a.^a,2._|_;)o, and a lim^^oo To(0 ^ (^^^ Appendix B), 
the two-point function {a^ax+i)i — > {a'^ax+i)^, U — >• 0. On the other hand, using the large 
U asymptotic formulae for limi,^oo to(0 (see Appendix B), we obtain: for [/ — > oo 

[Trail + l))-^sin7rpe(/ + 1) 
lim (a+a^+i)i PS <{ +(7ra(/ - l))-^sin7rpe(/ - 1), / = 2s + 1, s = 0, 1,..., (79) 
(1 + (-1)^) (7r0"^sin7rpe^, l = 2s, s = l,2,.... 



In particular, for pe = 1/2) hmi/^cx) hniL^oo(flx = 0. From now on, we consider only the 

case of Pe = 1/2. In the infinite-chain limit, limi^oo(a+aj;+2s)o = 0, s = 1, 2, . . ., thus 

, + > _ r limi,^oori(0, Z = 2s + 1, s = 0, 1, . . . , 

il'^^'^^^^+'^^-i (-l)^alim^_ooro(0, l-2s, s = l,2,.... ^^^^ 
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By ||2^ (see also Appendix B), 



and 



K f ( 

lim Tn[2s) = — / — — 



cos(2sA;) 



sin^ k 



-Ak 



2s+l 



(2.-1)! 



lim ri(2s + 1) 

L— »oo 



23S+2 g\ 



(2s- 1)!! 
23^+2 ^! 

2s + 1 



47r y Vl — /t^ cos^ k 



i^(s + ^,s + ^,2s + l;«' 



=dA; 



(81) 



f82) 



where F(a^b, c; z) stands for the hypergeometric function By means of an integral 

representation of F{a, b, c; z) p9|, to(2s) can be rewritten in the following form 



2s+l 

hm ro(2s) = (-1)^ -- 



exp 



s In 



1 - kH 



dt 



(83) 



which is suitable for determining the asymptotics of ro(2s) for s — > oo. Using the Laplace 
asymptotic formula ||30| for the integral in (p3D, we obtain 



lim To (2s) 

L— >oo 



1 



where the correlation length ^ is defined by 

^ = 



exp(-s/0 



1 



21n Jl 



C/\2 
4 / 



and has the following asymptotic behaviour: 



^ ~ for f/ ^ 0, 



and 



(85) 



(86) 



21n(f)' 



for f/ ^ oo. 



(87) 



In Fig. 1^ we show ^ as a function of U, obtained from the exact formula ([85| ) and from fitting 
numerically calculated correlation function {a^ax+2s)i with the formula const exp {—s / ^) . We 
plot also in Fig. the asymptotic formulae (|86D and (|87|). As a matter of fact, the formula 
( P5| ) approximates the exact value of ^ with accuracy 10"^^ up to f/ = 1. The formula ( [57| ) 
approximates ^ with accuracy 2, 5 ■ 10"^^ for U > 10. Consequently, we find the asymptotic 
forms of the two-point function (a+a^+;)i: 



I — '^s, I (a!^a^+2s)i| 



1/ 

4 



2tt 



1 + 



exp(-2s/0 
2s 



(88) 
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and 



/ = 2s 



57r 



exp(-l/0 



exp(-2s/0 
2s 



^9) 

Having obtained the correlation functions {n^)i and (a|^a^+;)i, we can easily calculate other 
correlations of interest. The short-range correlation functions (|2^) , (|23|) , for / = read: 

>^p(0) = Pe(l -Pe), 



5p(0) = --aro(0). 



Then, for / > and any pe < 1/2 we find 



C,{1) = (-l)'«Vo^(0) 



/ = 2s 

(a+a^+;)g + aV2(/), / = 2s, 



s = 0,l. 



1,2, 



(90) 
(91) 

(92) 



+ (-1) 



n 2 



- aT-o(O) 



/ = 2s + 1, 



0,1, 



s = l,2. 



(93) 



(a+a^+i)2 + aV2(/), / = 2s, 

In the infinite-chain limit and for pe = 1/2, the above short-range correlation functions assume 
the form: 



-lya' lim TolO) 

L— ►00 



limi^ooi"! (O5 



/ = 2s 



a^\imL^^T^{l), / = 2s, 



£(/) = lim Cp{l) 

-1, s = 0,l..., 
s = l,2.... 



(94) 



S{1) = lim Sp{l) 



-V. 



2s + 1, 
2s, 



0,1, 
1,2, 



(95) 



1 r ^n^l / limL_oorf(/), / 
--a2r^ro(0)J I 

The short-range correlations, given by ( P^ and (^), are shown in Fig. ^ as functions of the 
distance, and in Fig. ^ as functions of U. As U increases, the absolute value of C increases, while 
the absolute value of S decreases. At distance / = 1, a minimum, known as the correlation hole 



25|, is clearly visible. Then, at distances greater than a few lattice constants, the short-range 



correlations become periodic, due to the exponential decay of the two-point function {a^ax+i)i- 
Using (^), (0), and (^5]), we easily find the infinite-chain long-range correlation functions 
V{7t), C, and S. For pe = 1/2 they read: 



lim Vp{n) 

L—*oo 



a 



hm r2(0), 

_L— >oo 



(96) 



C = lim Cr 



-a' lim r2(0), 

4 L— +00 



S = ^(l-2a limro(O))' 
lb V L-+00 / 



(97) 



(98) 



For checkerboard configurations, 'P(vr) and C are proportional. Using the asymptotic formulae 
of Appendix B, we find that as f/ 00, V{7!') — > 1/4, £ ^ 1/16, and 5 — 0. The variation of 
long-range correlations versus U is shown in Fig. |^. As U increases, the quantum fluctuations 
decrease, but rather slowly, disappearing only in the limit of infinite U. For U = 5, C attains 
77% of its saturation value, and for U = 8.3, 90% of it. 
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5 Correlation functions — numerical results 



In the previous section we have been able to determine analytically the spatial and U de- 
pendence of a number of correlation functions, in the case of G = [1/2]. Here, our aim is 
to present the spatial and U dependence of the same kind of correlation functions, but with 
other choices of the set G. Specifically, we consider the following sets of low-period ground- 
state configurations: [1/6] i, [1/6] 2, [1/6] 3, [1/3], [2/5], [3/7], [1/4] 2; atomic, 2- molecular and 
3- molecular configurations are included in the list. The corresponding stability intervals of 
U, determined approximately by means of restricted ground-state phase diagrams obtained in 
T§, read: AU[ye]^ = [0.6, 00), AU^yeh = [0-3,2.7], AU^yeh = [0-l>l-0], At^[i/3]i = [0-2, 00) 



At^[2/5]i — [0-1) ^U[3/7]i — [0-05, 00), Af/[i/4]2 = [0.3, 2.2]. In contradistinction to the previ- 
ous section, the results reported here have been obtained numerically. We have used numerical 
procedures elaborated for a study of spin chains [0, with the main part being an exact diag- 
onahzation of the matrix {x\hb{w)\y) , defined by (^61) , (|47|) . The size of the chains has varied 
from a few hundred up to a few thousand sites. In the calculations, we have imposed the free 
boundary conditions on a chain whose number of sites, L, was chosen in such a way that L 
was a multiple of 4g, with q being the period of the configurations in G under consideration. 
The one-point, broken symmetry averages {nx)b,g, and the two-point, symmetric, short-range 
correlations, Cb,x{^) and Sb^xi^), have been calculated using the formulae (|55D. To weaken the 
effect of the chain ends, the site x and distance / have been chosen so that all the sites involved 
in the calculations of correlations were well inside the chain. The results obtained for chains of 
different sizes have been extrapolated to the infinite-chain limit. Practically, it has appeared 
that the results obtained for chains of L = 420 coincided with the corresponding infinite-chain 
results with accuracy (the relative error) exceeding 10~^^. 

In Figs. 1^, ^, we show short-range correlations as functions of distance. All the plots of 
two-point correlation functions, versus distance, exhibit a minimum at distance one, i.e. the 
so called correlation hole p5|, and after a few further steps become perfectly periodic, because 
of the exponential decay of the function (a'^ax+i) (see below). In Figs. |^ - |T2| we have plotted 
the short-range correlations, for a few small distances, as functions of U. Typically, in the 
presented cases the short-range correlations vary monotonically with U. However, S{1) for 
G = [1/6] 1 (shown in Fig. |T2[) exhibits a minimum for some U < 1. It is caused by a non- 
monotonic behaviour of {nx)g at nearest neighbours of the occupied site; for small U {U < 0.5 
roughly) it increases, and then decreases (see Fig. |T^). We note also that, as [/ — > 00, 5(0), 
(S(l) do not tend to 0, for molecular configurations. This is the case also for S{2) in the case 
of 3-molecular configurations [1/6] 3. The reason is that for molecular configurations, the local 
density of electrons at site x does not tend to sls U —>■ 00, but obeys, for molecules greater 
than 3-molecules, some nontrivial distribution. This is shown in Fig. |T^. 

Of particular interest is the spatial decay of the two-point correlation function \{a'^ax+i)b\'^- 
The exact calculations carried out for checkerboard configurations, in the previous section, 
suggest that for periodic ion configurations this decay is of the form /""^ exp(— //^), with some 
positive 7. Thus, this decay can be characterized by the correlation length ^ and the index 7. 
Specifically, for checkerboard configurations we have obtained ^ as a function of U (|5^, with 
= U/2, for small U, that is amounts to a half of the gap at the Fermi level, and with 
7 = 1, according to (pSj), (pOf). For G other than the checkerboard configurations, we have 
approximated the large distance behaviour of \ {a^ax+i)b\'^ by the formula const exp{—l / ^) , thus 
ignoring the possible power-law prefactor The data are shown in Fig. . Comparing the 
data obtained in such calculations of ^ with the exact results in the case of the checkerboard 
configurations, we estimate the error to be at the level of a few per cent. To reach such an 
accuracy, in the range of U < 0.4 it was necessary to consider a chain with L = 2 ■ 10^, while 
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for larger U it was enough to take L = 10^. For all the studied periodic configurations, we have 
found that for small U the inverse correlation length ^ constA, with A being the gap at 
the Fermi level, and A ~ f/ . 

Finally we present our results concerned with the defined in Section 3 long-range correlation 
functions. Our purpose was not only to calculate the correlations Vfi^k), Cf, and Sj, in the 
infinite-chain limit, but also to observe their dependence on the chain size L, in a finite chain 
with the free boundary conditions imposed {b = /). To calculate the long-range correlations, 
we used their definitions (0), (|3^) , (^31), which involve all sites of the chain, and the simplified 
expressions, involving L/2 central sites of the chain, 

^K^) = ill ^""'^lAl), (99) 

1=0 



2 ^^^^^ 

^'f=L^ (100) 

1=0 



2 



1=0 



where x = L/A. As shown in Fig. |15[ for particular sets G and values of U, the long-range 
correlations vary linearly with 1/L, and the simplified expressions, being less sensitive to the 
boundaries of the chain, for given L constitute a lot better approximation to the infinite-chain 
value of correlations than the defining formulae. The dependence of infinite-chain values of 



V, C, and S on U, for different choices of the set G, is displayed in Figs. 0, [T^. As a rule, 
£ increases towards a positive saturation value, while S decreases towards zero, for all the 
considered sets G, except G = [l/6]3, when S exhibits a minimum at f/ ^2. For U > 2, 

5 increases towards a positive asymptote. An inspection of the local density of electrons, as 
a function U, reveals a non-monotonic behaviour of this density at nearest neighbours of the 
central occupied site of the 3-molecule: it has the maximum at f/ ^ 2 and decreases towards 
a positive asymptote, a.s U —>■ oo (see Fig. 0). Therefore, the non-monotonic behaviour of the 
local density is responsible for the non-monotonic behaviour of S. 

6 Summary and discussion 

In the paper, we have studied analytically and numerically, the spatial and U dependence 
of short-range and long-range ground-state correlation functions of electrons in the Falicov- 
Kimball chain. We have chosen the electron- ion coupling U > 0, which amounts, in our 
case, to attraction between the electrons and ions. The obtained results can be transformed 
to the case of repulsion by means of hole-particle transformations P|. We emphasize also 
that only nearest-neighbour, uncorrelated with ions, hopping of electrons has been taken into 
account. The properties of electron correlations in the case of extended or correlated with ions 
hopping will be presented elsewhere. The calculations have been carried out in different ground 
states, specified by the sets G of periodic ground-state ion configurations and the corresponding 
stability intervals AUg- Typically, for a specified G we considered an interval of U having only 
a non-void intersection with the stability interval of G, AUg- 
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Outside the stability intervals, the results obtained refer to tight-binding electrons in peri- 
odic external potentials, given by G. 

We mention that, due to the Jordan- Wigner transformation, our results can also be in- 
terpreted in terms of spin-1/2 isotropic XY chains, in periodic transverse external magnetic 
fields. 

It is well known that the two-point correlation function {a'^ax+i)b, / = 0, 1, . . ., amounts to 
the one-body reduced density matrix in the ground state. This quantity enables one to calculate 
ground-state averages of one-body observables and is therefore of crucial importance for 
all methods that describe many-body systems by means of single-particle formalism, like the 
density functional theory, for instance. The spatial decay of the one-body reduced density 
matrix determines the degree of locality of many quantities, which are relevant for describing 
properties of metals and insulators. For this reason, its asymptotic long-distance behaviour has 
been the subject of interest for long time. Quite recently, the spatial decay of the one-body 
reduced density matrix was investigated in some tight-binding models of metals and insulators 
(see |3^, and the references quoted there). In the authors considered a model 

of an insulator with two-bands separated by a gap, due to two bare electronic states, with 
different bare energies, at each site, with equal in-band transfer integrals and a weak interband 
hybridization. The asymptotic long-distance behaviour of the single-particle density matrix has 
been studied, in dimensions D = 1,2, and 3. The dependence of the correlation (decay) length 
on the energy parameters of the model and the existence of the power-law prefactor, of the form 
l-D/2^ has been demonstrated. The results obtained in our paper can be viewed as the results 
referring to models of many-band insulators, where the bands separated by gaps arise due to 
a periodic spatial modulation of the bare energies. In the case of two-bands (Section 4), we 
have been able to calculate exactly the correlation functions, and then extract, in particular, 
the asymptotic long-distance behaviour of the single-particle density matrix. In the considered 
here case oi D = 1, we have found the prefactor, as in [Q, and the exact expression 
for the correlation length ^, given by (P5|), as a function of the periodic-potential strength U 
(which amounts also to the gap width). For a weak potential, the inverse correlation length 
f» f//2, while for a strong potential, ^ 21n([//2). 

In ||3^, Macedo et al calculated the correlation function S{1), given by (p5|), as the T 
limit of the corresponding grand-canonical correlation function, that is, with the ground-state 
symmetric average in definition (|23| ) replaced by the grand canonical average. The chemical 
potential was chosen in such a way that the ground-state was characterized by the set [1/2] i 
of checkerboard ion configurations. In the case of the one- dimensional system under consider- 
ation, the grand canonical correlation functions should coincide with the canonical ones. The 
results in have been obtained using the method of small-cluster exact diagonalization and 
extrapolation techniques to the infinite chain, where the cluster size was limited to at most 10 
sites. To compare the values of S{1), obtained from the exact formula (P^D, and those displayed 



in figures of one has to change the sign of U in (p5D, and to multiply our S{1) by 3/4. The 



results of exact calculations are shown in Table 1, and they differ qualitatively from the data 
displayed in p5|. 

The ground-state phase diagram of the Falicov-Kimball model is reach in quantum phase 
transitions, that is the phase transitions where the role of thermal fiuctuations is played by 
quantum fiuctuations ||3^. These transitions are driven by such control parameters as pe, Pi, 
and U. They occur, when for some critical value of a control parameter, the nature of the ground 
state changes. A detailed description of quantum phase transitions in a many-body, interacting 
system, like the Falicov-Kimball model, is a hard task (for recent results see |^, |^). We 
would like to point out that the analytic results of Section 4, can be used to describe a simple 
instance of a quantum phase transition, driven by U, with the critical value Uc = 0. Consider an 



18 



Table 1: 35(/)/4 



/ 


U = 


f/ = -4 


U = -8 





0.375 


0.688 


0.729 


1 


-0.236 


-0.658 


-0.719 


2 


0.188 


0.630 


0.710 


3 


-0.196 


-0.632 


-0.710 



infinite chain, whose particle densities are set to 1/2, pj = pe = 1/2. Obviously, as U = Uc = 
there is no order in the ion subsystem, whereas the electrons are distributed uniformly. There is 
no gap in the electron energy spectrum and the truncated electron density- density correlation 
function, {rixUx+i) — 1 / 4:, exhibits a power-law decay (with the oscillations owing to sin2(7r//2) 
imposed) as / — cxd, according to (^). On the other hand, for an arbitrary nonzero U, the ion 
subsystem becomes checkerboard ordered @, |^ . The electrons follow the periodic distribution 
of the ions, for instance the broken-symmetry average {nx)i is modulated with period 2 (|67D. 
A gap A = U appears at the Fermi level, and correspondingly, the large-distance behaviour 
of the broken-average truncated density-density correlation function, {n^n^+i)! —{nx)i{nx+i)i, 
changes to an exponential one, of the form /~^exp(— //^) (see (|8]), (|89|) ), with the correlation 
length ^ ~ U^^, which diverges as f/ ^ 0. 
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7 Appendix A 

Here we define the short- and long-range correlation functions for the ion subsystem. For any 
ion configuration g & G, and its restriction to a finite chain whose number of sites, L, is a 
multiple of the period q of the ground-state configurations g, 

L-1 

q-'Y.9- = P^ = L-'Y.9- (102) 

geG x=0 

The short-range ion-ion correlation function is 

Eb,x{l) = (^^- ~ (^^+' ~ ' (10^) 

which for b = p is translation invariant. The degree of order can be "measured" either by the 
ion-order parameter, If, g{k), which for given g & G and k = 27rpe = 2Tip/q is 

L-1 

h^g{k) = 5^e*^- {g, - p,) ; ^^/^(fc) = 0, (104) 

x=0 g 

or by the long-range correlation, '^geG \^b,g{k)\'^, 

L-1 

E I^M(^)r = E e^^(^-^)g-i E [{g^ - p,) {g, - p,)] . (105) 

<?eG x,y=0 geG 



19 



The function k ^ q ^ ^^^q \Ib,g{k)\'^ is the static structure factor for the ions. For periodic 
boundary conditions, 

q-l 

q-' J2 Up,sik)\' = = q-' Yl ^^'^''"'^ - P^) iOy ' P^) > 0, (106) 

g&G X, y=0 

where the above inequahty follows from the argument presented in the Appendix of |^ . Thus, 
the square of the absolute value of the ion-order parameter is a measure of the ion-ion long-range 
order. 



8 Appendix B 

For pe < 1/2) the large U asymptotic behaviour of limi^ooTo(0 is smooth in and reads: 

f TT r ^ f sin(7rpeO 

for {/ — > cxD, lim ro(/)^ — < ; 

L^oo na I I 

sin(7rpe0 sin(7rpe(^ + 2)) sin(7rpe(Z — 2)) 



2a2 



1 + 2 



1-2 



2a^ 



sin(7rpe/) _ 5 / sin(7rpe(/ + 2)) sin(7rpe(/ - 2)) 
/ 2 V / + 2 1-2 

^ sin(7rpe(/ + 4)) ^ sin(7rpe(/ - 4)) 



/ + 4 



I- A 



In particular. 



for [/ ^ oo, 



lim ro(0)^Pe/a. 



L^oo 



(107) 



(108) 



On the other hand, the small U asymptotic behaviour of lim2,^ooTo(0 is singular at pe = 1/2. 
For Pe < 1/2, the following recursive formula holds: 



for f/ ^ 0, 
Thus, knowing that 

for ^ 0, 



,,, sin(7rpe(^ — 1)) 
lim ro(0^ ^^.7^ " - lim r,{l - 2) 



lim ro(0)^— In | tan [7r(l + 2pe) /4] 

L^oo 2tX 



(109) 



;iio) 



we can obtain the asymptotic formula for any even I > 2: 

i-ir 



lim ro(2n) 

L— >oo ZTT 



+ — 



2j-l 



n-l 



4^(-l)^(n-j)(^-J + l) 



sin(2j - l)7rpe 



2j-l 



+ {2n' - 1/2) In tan [7r(l + 2pe)/4] - - 



1 sinvrpe 



2 COS^ 71 Pe 



+ ... 



Then, for pe = 1/2 and / = 0, 2 we find: 
for f/ ^ 0, 



lim ro(0)^ — (ao + 21n2 - Ina), 

L^oo 27T 



111) 



'112) 
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for U^O, lim ro(2)^— (2-ao-21n2 + lna), (113) 

L^oo 2,TT 



with ao given by 

1 

fin = — 

2^ n(2n)!! 

n=l 



ao = -y%^- 0.687. (114) 



9 Appendix C 

The relation between the integral in (^) and the hypergeometric function can be obtained as 
follows: 



K /" ( 

iT^"°(2.) = -y-^ 



cos(2sfc) 

ak 





oo 



.2 sin^ k 



— [dk COs{2sk) y i^:L^^2i .^2j ^ 





oo 



= ILy /■ cos(2sA;) sin^J' A; 

4vrf^ (2j)!! J ^ ' 

_ . ^ ,, (2j-l)!! vr (2j)! 

4vrf^ (2j)!! ^ ^ 22.- (j + .)!(j - .)! 

^ K ^ (2j - 1)!! (2j)! 

J — * 

OO 



[(2j - 1)!!]^ /t^J' 
(j + .)!(j-s)! W 

[(2s + 2j - 1)!!]2 1 



4 ~^ (2s + j)' j' 22(«+J) 

= f-l?^V^^^^±i±^ — 
^ ^ 4 A.r(2s + l+j) j! 

^ ' 4n r(2s + l) ^ 2' 2' ' ^ 

, ,,k2^+^ (2s- 1)!! ,11 2, , , 



The Laplace asymptotic formula for integrals reads 



1 



2n 



where 



/(t)e^-Wdt ™ \/-;5%y/(^o)e^^(*»), (116) 



5(to) = max5(t). (117) 
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In the case of the integral in the last line of (|8 



1 ^ , K^-t) 

, , =ln-^ ^, 

,/t{l - t){l - kH) 1-kH 



with 

1 - VT^^ 



(119) 



and 



S{to) = \n{4) = In 



0) 



(1 - VT^)VT^' 



S"ito) = ^""^ (120) 



Consequently, the asymptotic form of the considered integral is 



TT 1 



i-VT 



2s 



.V^^e-^/^ (121) 



with ^ given by 
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Figure 1: Continuous lines: limi_+oo ^olO' given by (|69D , versus U, for Z = 0, 2 and different 
values of pe- Dashed lines: the small U asymptotes. Dotted lines: the large U asymptotes. 
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Figure 2: Continuous line: the exact correlation length for the checkerboard configurations, 
given by (|^), versus U. Dashed line: the small U asymptote (^61). Dotted line: the large 
U asymptote (W^ - Open circles: ^ obtained from fitting numerically calculated correlation 
function {a'^ax+2s)i with the formula const exp(—s/^). 
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Figure 3: (a) C{1), given by (0), versus /, for G = [l/2]i, and different values of U. (b) 
given by (|95D , versus /, for G = [l/2]i, and different values of U. Filled triangles with base up: 
U = 5, filled triangles with base down: U = 2, filled circles: U = 1, filled squares: U = 0.5. 
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Figure 4: (a) given by (@), versus U, for G = [l/2]i, and different values of /. (b) S{1), 
given by (P^D, versus [/, for G = [l/2]i, and different values of /. In the scale of the figure, the 
plots for odd / > 5 cannot be distinguished from the plot for / = 3. Similarly, the plots for 
I > 4 cannot be distinguished from the plot for I = 2. 
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Figure 5: Continuous line: £ for G = [1/2] i, given by (p7|), versus U. Dashed line: S for 
G = [1/2] 1, given by (]98|), versus U. 
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Figure 6: Infinite chain limit of Cf^x{l), given by (^2]), versus /, for G = [l/3]i (a), G = [l/4]2 
(b), and for different values of U. Infinite chain limit of S{1), given by (p3D, versus /, for 
G = [l/3]i (c), G = [l/4]2 (d), and different values of U. Filled triangles with base up: U = 5, 
filled triangles with base down: U = 2, filled circles: U = 1, filled squares: U = 0.5. 
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Figure 7: Infinite chain limit of Cf^x{l), given by (^2]), versus /, for G = [2/5]i (a), G = [3/7]i 
(b), and for different values of U. Infinite chain limit of S{1), given by (p3D, versus /, for 
G = [2/5]i (c), G = [3/7]i (d), and different values of U. Filled triangles with base up: U = 5, 
filled triangles with base down: U = 2, filled circles: U = 1, filled squares: U = 0.5. 
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Figure 8: Infinite chain limit of Cf^^i^), given by (p2D, versus Z, for G = [l/6]i (a), G = [l/6]2 
(b), and for different values of U. Infinite chain limit of S{1), given by (pSf), versus I, for 
G = [1/6] 1 (c), G = [l/6]2 (d), and different values of U. Filled triangles with base up: U = 5, 
filled triangles with base down: U = 2, filled circles: U = 1, filled squares: U = 0.5. 
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Figure 9: Infinite chain limit of Cf^rc{l), given by (^), versus U, for different distances I, and 
different sets G. 
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Figure 10: Infinite chain limit of Sf^^i^), given by (^), versus U, for different distances /, and 
different sets G. Whenever a plot for G = [3/7] i is missing, it is indistinguishable from the plot 
for G = [2/5] 1, in the scale of the figure. 
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Figure 11: Infinite cliain limit of Cf^^i^), given by (p2|), versus U, for different distances /, and 
different sets G. Wfienever a plot for G = [l/Qjs is missing, it is indistinguisliable from tlie plot 
for G = [l/6]2, in the scale of the figure. 
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Figure 12: Infinite chain limit of S{1), given by (P^D, versus U, for different distances /, and 
different sets G. 
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Figure 13: The local electron density in a fixed ion configuration g, versus position x, as a 
function of U. The ion configuration g = 100000 for G = [l/6]i, c/ = 110000 for G ^ [1/6] 2, 
and g = 111000 for G = [1/6)3. 
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Figure 14: Numerically calculated correlation length ^ versus U, and numerically calculated 
const A~^, with const adjusted so that ^ pa const A~^, for different sets G. The ^ has been 
obtained from approximating the large-distance behaviour of \{a'^ax+i)b\'^ by constexp{—l/^). 
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Figure 15: Vf(27ipe), given by (|3T|), and V'j-{27ipe), given by (^9]), versus L ^, for U = 2, and 
G = [l/2]i (a), G = [3/7]i (b). Filled squares: P/(27rpe), filled circles: V'f{2np^). 
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Figure 16: (a) P(27rpe), given by (|43|) , versus ?7, for different sets G, (b) £, given by 
versus [/, for different sets G, (c) S, given by (0), versus [/, for different sets G. In (a) and 
(c) tlie plot corresponding to G = [3/7]i is missing, since it is indistinguisliable from the plot 
for G = [2/5] 1. 
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Figure 17: (a) V{2Txpe), given by (|43|) , versus t/, for different sets G, (b) £, given by 
versus [/, for different sets G, (c) iS, given by (|^), versus [/, for different sets G. 
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